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Abstract 

We demonstrate a model-independent method of estimating qualitative dynamics of an accel- 
erating universe from observations of distant type la supernovae. Our method is based on the 
luminosity-distance function, optimized to fit observed distances of supernovae, and the Hamilto- 
nian representation of dynamics for the quintessential universe with a general form of equation of 
state p = w(a(z))p. Because of the Hamiltonian structure of FRW dynamics with the equation 
of state p = w(a(z))p, the dynamics is uniquelly determined by the potential function V(a) of 
the system. The effectiveness of this method in discrimination of model parameters of Cardassian 
evolution scenario is also given. Our main result is the following, restricting to the flat model 
with the current value of Q m ,o = 0.3, the constraints at 2a confidence level to the presence of p n 
modification of the FRW models are —0.50 < n < 0.36. 
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I. INTRODUCTION 

Recent observations of type la supernovae |3J] supported by WMAP measurements of 
anisotropy of the angular temperature fluctuations 4| indicate that our Universe is spatially 
flat and accelerating. On the other hand, the power spectrum of galaxy clustering 5( 
indicates that about 30% of critical density of the universe should be in the form of non- 
relativistic matter (cold dark matter and baryons). The remaining, almost two thirds of 
the critical energy, may be in the form of a component having negative pressure (dark 
energy). Although the nature of dark energy is unknown, the positive cosmological constant 
term seems to be a serious candidate for the description of dark energy. In this case the 
cosmological constant A and energy density e\ = A/(8ttG) remain constant with time and 
the corresponding mass density p\ = £\/c 2 = 6.44 x 10~ 30 (Q\/0.7)(h/0.7)g cm~ 3 , where h is 
the Hubble constant H expressed in units of 100 km s _1 Mpc _1 and = 0.7±0.1. Although 
the cold dark matter (CDM) model with the cosmological constant and dust provides an 
excellent explanation of the SNIa data, the present value of A is ~ 10 123 times smaller than 
value predicted by the particle physics model. Many alternative condidates for dark energy 
have been advanced and some of them are in good agreement with the current observational 
constraints . Moreover, it is a natural suggestion that A-term has a 

dynamical nature like in the inflationary scenario. Therefore, it is reasonable to consider 
the next simplest form of dark energy alternative to the cosmological constant (w = —1) 
for which the equation of state depends upon time in such a way that p = w(a(z))p, where 
w = p/p is a coefficient of the equation of state parametrized by the scale factor or redshift. 
It has been demonstrated 3, ^| that dynamics of such a system can be represented by 
one-dimensional Hamiltonian flow 

W = y + V(a), (1) 

where the overdot means differentiation with respect to the cosmological time t and V(a) is 
a potential function of the scale factor a given by 

V{a) = -*^a a , (2) 
o 

where p c fr is the effective energy density which satisfies the conservation condition 

a . 

PcS = - 3 - (p cff +PeS)- 
CL 



For example, for the ACDM model we have 



p cfr = A + p mfl a 
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PcS = -A + 0. 



(3) 



Of course the trajectories of the system lie on the zero energy surface Ti, = 0. 

Hamiltonian can be rewritten in the following form convenient for our current recon- 
struction of the equation of state from the potential function V, namelly 



where p x = dL/dx, here overdot means differentiation with respect to some new 
reparametrized time t — > r: dr = \H \dt, x = a/a = 1 + 2. 

For example, for mixture of non- interacting fluids potential V(x) takes the form 



where pi = WiPi for z-th fluid and Vz Wi = const, (similar to the quiessence model of dark 
energy) . 

Due to the Hamiltonian structure of Friedmann-Robertson- Walker dynamics, with the 
general form of the equation of state p = w(a(z))p, the dynamics is uniquely determined 
by the potential function V(a) (or V(x)) of the system. Only for simplicity of presentation 
we assume that the universe is spatially flat (in the opposite case trajectories of the system 
should be considered on the energy level Ti = \Qk,o)- 

Let us note that from the potential function we can obtain the equation of state coefficient 
w = p/p, 



The term d(ln V)/d(hia) has a simple interpretation as an elasticity coefficient of the po- 
tential function with respect to the scale factor. 

Thus from the potential function V both p(a) and p can be unambiguously calculated 



n(p,x) = ^ + v(x), 



(4) 




(5) 




(6) 



p(a) 
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V(a) 




(7) 
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II. THE EXPANSION SCENARIO FROM SUPERNOVAE DISTANCES 



As it is well known in a flat FRW cosmology the luminosity distance Dl and the coordi- 
nate distance r to an object at redshift z are simply related as 

t dt ' (rt 

(c = 8nG = 1 here and elsewhere). From equation (JBJ) the Hubble parameter is given by 

d (DAz 



H(z) = [In a] 



dz 



(9) 



It is crucial that formula (jUJ) is purely kinematic and depends neither upon a microscopic 
model of matter, including the A-term, nor on a dynamical theory of gravity. Due to 
existance of such a relation it would be possible to calculate the potential function which is: 



Via) = --H 2 a 2 = - ldz n )^ +z [i . (10) 



1 tt2 2 Idz \ l+z /J 
— tl (1 = ; r 

2 2(1 + z) 

This in turn allows us to reconstruct the potential V(z) from SNIa data. Let us note 
that V(z) depends on the first derivative with respect to z whereas w(z) is associated with 
the second derivative. 

Let us also note that from of the potential function for a one-dimensional particle- universe 
moving in the configurational a (or x)-space can be reconstructed from recent measurements 
of angular size of high-z compact radio sources compiled by Gurvits et al. jlfij ]. The corre- 
sponding formula is 



v(a) = z:\ :,: ,n , (n) 



[£( D A (z)(l + z))\ 
2(1 + z) 

where the luminosity distance Dl and the angular distance Da are related by the simple 
formula 

D L {z) = {l + zfD A {z). (12) 

Since the potential function is related to the luminosity function by relation (jlOj) one 
can determine both the class of trajectories in the phase plane (a, d) and the Hamiltonian 
form as well as reconstruct the quintessence parameter w(z) provided that the luminosity 
function Dl(z) is known from observations. 

Now we can reconstruct the form of the potential function (HUB using a natural Ansatz 
introduced by Sahni et al. |16j. In this approach dark energy density which coincides with 
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p c s is given as a truncated Taylor series with respect to u — (1 + z) 

PDE = PeflF = A + Aiu + A 2 u 2 . 



This leads to 



H(z) 



H oy /n mfl (l + zf + A + + 2) + A 2 (l + 



(13) 



and 




^o7o v/fi m , (l + zf + A + A 1 (l + z)+ A 2 (l + ^) 2 ' 




(14) 



The values of three parameters Ao,Ai,A 2 can be obtained by applying a standard fitting 
procedure to SNIa observational data based on the maximum likelihood method. 
The potential function (fTU|) written in terms of (1 + z) is 



where Ai/(3H%) = Q i>0 , (i = A,t, k). 

Our approach to the reconstruction of dynamics of the model is different from the stan- 
dard approach in which w(z) is determined directly from the luminosity distance formula. 
It should be stressed out that the latter approach has an inevitable limitation because the 
luminosity distance dependence on w(z) is obtained through a multiple-integral relation that 
loses detailed information on w(z) |l2|. In our approach the reconstruction is simpler and 
more information on w(z) survives (only a single integral is required). Our approach is also 
different from the concept of reconstruction of potential of scalar fields considered in the 
context of quintessence |8( . 

The key steps of our method are the following: 

1) we reconstruct the potential function V(a) for the Hamiltonian dynamics of the 
quintessential universe from the luminosity distances of supernovas type la; 

2) we draw the best fit curves and confidence levels regions obtained from the statistical 
analysis of SNIa data; 

3) we set the theoretically predicted forms of the potential functions on the confidence 



V(a(z)) 





or in dimensionless form 




(16) 
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levels diagram; 

4) those theoretical potential which escape from the 2a confidence level is treated as being 
unfitted to observations; 

5) we choose this potential function which lie near the best fit curve. 



Our reconstruction is an effective statistical technique which can be used to compare a 
large number of theoretical models with observations. Instead of estimating some revelant 
parameters for each model separately, we choose a model-independent fitting function and 
perform a maximum likelihood parameter estimation for it. The obtained confidence levels 
can be used to discriminate between the considered models. In this paper this technique is 
used to find the fitting function for the luminosity distance. 

The additional argument which is important when considering the potential V(a) is that 
it allows to find some modification in the Friedmann equations along the "Cardassian expan- 
sion scenario" Q|. This proposition is very intriguing because of additional terms, which 
automatically cause the acceleration of the universe 0,0, HQ- These modifications 
come from the fundamental physics and these terms can be tested using astronomical ob- 
servation of distant type la supernovae. For this aim the recent measurements of angular 
size of high-redshift compact radio sources can also be used Q| . 

The important question is the reliable data available. We expect that supernovae data 
would improve greatly over next few years. The ongoing SNAP mission should gives us 
about 2000 type la supernovae cases each year. This satellite mission and the next planned 
ones will increase the accuracy of data compared to data from the 90s. In our analysis we 
use the availale data starting from the three Perlmutter samples (Sample A is the complete 
sample of 60 supernovae, but in the analysis it is also used sample B and C in which 4 and 
6 outliers were excluded, respectively). The fit for the sample C is more robust and this 
sample was accepted as the base of our consideration. For technical details of the metod the 



reader is referred to our previous two papers 
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In Fig. 1 we show the reconstructed potential function obtained using the fitting values 
of Ai as well as A4. The red line represents the potential function for the best fit values of 
parameters (see TABLES ITU IIHI and IIV |) . In each case the coloured areas cover the confi- 
dence levels 68.2% (la) and 95.4% (2a) for the potential function. The different forms of 
the potential function which are obtained from the theory are presented in the confidence 
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levels. Here we consider one case, namely the Cardassian model. In this case the standard 
FRW equation is modified by the presence of an additional p n term, where p is the energy 
density of matter and radiation. For simplicity we assume that density parameter for radi- 
ation is zero (see TABLE. QJ). The Cardassian scenario is proposed as an alternative to the 
cosmological constant in explaining the acceleration of the Universe. In this scenario the 
the Universe automaticaly accelerates without any dark energy component. The additional 
term in the Friedmann equation arises from exotic physics of the early universe (i.e., in the 
brane cosmology with Randall-Sundrum version n — 2). 



TABLE I: The forms of the potential functions in dimensionless form for two cases: ACDM model 
and Cardassian scenario. 





The form of the potential 
function 


Position of the maximum 
predicted by theorethical 
models 


Position of the maximum 
from reconstruction 










sample 


value 










A 


0.613 


Perlmutter model 


v(x) = -i^oi" 1 + n A . x 2 ] 


3/1 n ™,o 
x o y 2 i-n m , 




C 


0.594 










A 


0.598 


Cardassian scenario 


V(x) = -^[Qm.,03;- 1 + n Cfi x 2 - 3 "] 




3(n-l) 


C 


0.599 



The solid lines in Fig. represent the potential functions for the different Cardassian 
scenarios of an accelerated expansion of the universe. 

For visualization the quality of performed fitting procedure the Hubble diagrams are 
presented on Fig. [21 

The general conclusion from our statistical analysis is that the potential function has 
maximum at some z = z (or a = a ) and that 

d 2 V 

a? w < °- (17) 

In the next section it will be demonstrated that the reconstructed potential function 
contain all informations which we need for reconstruction of the dynamics of any cosmological 
model on the phase plane. 
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FIG. 1: Confidence levels for the potential function V for the model with a) ^3 = 0.3, (sample 
C), b) A 3 = 0.05, (sample C), c) A 3 = 0.3, (sample A), d) A 3 = 0.05, (sample A); and for the 
Perlmutter model e) (sample C), f) (sample A). The best fit for each model is represented by the 
red curve. Black and blue curves represent potential functions for theoretical Cardassian models 
for n mfi = 0.3 (a,b,e,f), Q m , = 0.05 (c,d) and n = 1/2, 1/3, 0, -1/3, -2/3, -4/3. 
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a) 



b) 




c) 



SNIadata h 
EdeS model 

Perlmutter model - 

O m = 0.3 fixed - 

n m(1 = 0.05 fixed - 

Lf,o'=0-7.n = 1/3 - 

n rn = 0.7, n = - 



0.2 0.3 0.4 0.5 




SNIa data 
EdeS model 
Perlmutter model 

n mQ = 0.3 fixed 
n m0 = 0.05 fixed 

fl G '= 0.7, n= 1/3 

ri r: = 0.7. n = 




d) 




FIG. 2: Hubble diagrams for Perlmutter SNIa data for sample A (a,b) and C (c,d). Coloured curves 
represent the best fitted magnitude-redshift relations (RGB colour standard) (The Perlmutter 
model and polynomial fits with fixed £l m ,o parameters respectively) and hypothetical Cardassian 
scenarios (CMYK colour standard) for different values of parameters £lc,o an d n. 

III. APPLICATION TOWARDS THE RECONSTRUCTION OF PHASE POR- 
TRAIT AND MODEL PARAMETERS 



The dynamics of the considered cosmological models is governed by the dynamical system 



dx 

x = — = y, 
dt y 



dy 
V =dt 



dV 

dx ' 



(18) 



with the first integral for (|18|) 7i = <^ x 2 /2 + V(x) = 0. The main aim of the dynamical 
system theory is the investigation of the space of all solutions (JTHj) for all possible initial 
conditions, i.e. phase space M. In the context of quintessential models with the equation of 
state p = w(a(z))p there exists a systematic method of reducing Einstein's field equations 
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TABLE II: Fitting results of the statistical analysis from distant type la supernovae data. The 
best fit values of parameters A{ and A4 for the polynomial fit (|14|) for the sample A of Perlmutter 
SNIa data. In the table the values of the best fit parameters as well as the most probable values 
are also given. The fit for the sample C is more robust one, because the two very likely reddened 
supernovae SN1996cg and SN1996cn have been removed 3] 





A 


A x 


A 2 


A 3 


M 


x 2 


Zmax 


Vmax 


best fit 


0.61 


0.09 


0.00 


0.30 


-3.38 


96.03 


0.659 


-1633.21 


max(P) 


u.ou_ 043 


09+ - 43 


0.00 


0.30 


O OQ+0.06 

°- oo -0.06 








best fit 


-0.25 


1.20 


0.00 


0.05 


-3.37 


96.84 


3.677 


-1011.88 


max(P) 


— 0.26±g;|| 


1 20+ - 44 


0.00 


0.05 


O q fi + 0.06 

°- ou -0.06 









TABLE III: Fitting results of the statistical analysis from distant type la supernovae data. The 
best fit values of parameters A{ and Ai for the polynomial fit (|14[) for the sample C of Perlmutter 
SNIa data. 





A 


Ai 


A 2 


A 3 


M 


x 2 


Zmax 


V m ax 


best fit 


0.74 


-0.04 


0.00 


0.30 


-3.43 


53.28 


0.676 


-1567.33 


max(P) 


7S+ - 40 


04+ - 44 


0.00 


0.30 


q ^q + 0.06 

°-^°-0.06 








best fit 


-0.13 


1.08 


0.00 


0.05 


-3.42 


53.83 


3.522 


-968.74 


max(P) 


-0 14+ - 40 


l.U8_o_ 39 


0.00 


0.05 


q 4 -,+0.06 









to the form of the dynamical system (fTHjl |l3[ . One of the features of such a representation 
of dynamics is the possibility of resolving of some cosmological problems like the horizon 
and flatness problems in terms of the potential function V(x). 

The phase space M (or state space) is a natural visualization of the dynamics of any 
model. Every point P = (x,y) G M corresponds to a possible state of the system. The 
r.h.s of the system ()18|) define a vector field F(P) = [y, — f^-] belonging to the tangent space 
TpM. Integral curves of this vector field define one-parameter group of diffeomorphisms 
0(P) called the phase flow. In the phase space the phase curves (orbits of the group 4>{P)) 
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TABLE IV: Fitting results of the statistical analysis for the Perlmutter model and both (A and C) 
samples of Perlmutter SNIa data 





sample 


^A,0 


^m,0 


M 


x 2 


Zmax 


Vmax 


best fit 


C 


0.71 


0.29 


-3.43 


53.29 


0.684 


-1573.69 


max(P) 


C 




3n+ 05 

u ' ou -0.05 


-3.43 








best fit 


A 


0.69 


0.31 


-3.39 


96.00 


0.632 


-1628.31 


max(P) 


A 


0-681851 


n oq+0.05 
u - oz -0.05 


-3.39 









represent the evolution of the system whereas the critical points y — 0, ^ = are singular 
solutions - equilibria from the physical point of view. The phase curves together with critical 
points constitute the phase portrait of the system. 

Now we can define the equivalence relation between two phase portraits (or two vector 
fields) by the topological equivalence, namely two phase portraits are equivalent if there 
exists an orientation preserving homeomorphism transforming integral curves of both sys- 
tems into each other. Following the Hartman-Grobman theorem, near hyperbolic critical 
points (V i ReAj 7^ 0, where Aj is the appropriate eigenvalue of linearization matrix A of the 
dynamical system) is equivalent to its linear part 



x = y, 

' d 2 V 



dx 2 

In our case the linearization matrix takes the form 



(x-x ). (19) 

(x=xq,0) 



A 



0. 



dx 2 ' u 



(20) 

J (x ,0) 

Classification of critical points is given in terms of eigenvalues of the linearization matrix 
since the eigenvalues can be determined from the characteristic equation A 2 — (Tr*4)A + 
detA = 0. In our case Tr^l = and eigenvalues are either real if fj^rL,^ ^ < or purely 
imaginary and mutually conjugated if fj^L^Q) > 0. In the former case the critical points 
are saddles and in the latter case they are centres. 
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The advantage of representing dynamics in terms of Hamiltonian (JTJ is the possibility to 
discuss the stability of critical points which is based only on the convexity of the potential 
function. In our case the only possible critical points in a finite donain of phase space are 
centres or saddles. 

The dynamical system is said to be structurally stable if all other dynamical systems 
(close to it in a metric sense) are equivalent to it. Two-dimensional dynamical systems on 
compact manifolds form an open and dense subsets in the space of all dynamical systems 
on the plane 2^|. Structurally stable critical points on the plane are saddles, nodes and 
limit cycles whereas centres are structurally unstable. There is a widespread opinion among 
scientists that each physically realistic models of the universe should possess some kind 
of structural stability - because the existence of many drastically different mathematical 
models, all in agreement with observations, would be fatal for the empirical method of 



science 
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Basing on the reconstructed potential function one can conclude that: 

1) since the diagram of the potential function V(a(z)) is convex up and has a maximum 
which corresponds to a single critical point - the quantity ^JrL^ ^ < at the critical point 
(saddle point) and the eigenvalues of the linearization matrix at this point are real with 
oposite signs; 

2) the model is structurally stable, i.e., small perturbation of it do not change the structure 
of the trajectories in the phase plane; 

3) since a = — ^ one can easily conclude from the geometry of the potential function that 
in the accelerating region (a > 0) V(a) is a decreasing function of its argument; 

4) the reconstructed phase portrait for the system is equivalent to the portrait of the model 
with matter and the cosmological constant. 

By an inspection of the phase portraits we can distinguish four characteristic regions in 
the phase space. The boundaries of region I are formed by a separatrix coming out from 
the saddle point and going to the singularity and another separatrix coming out of the 
singularity and approaching the saddle. This region is covered by trajectories of closed and 
recolapsing models with initial and final singularity. 

The trajectories moving in region IV are also confined by a separatrix and they correspond 
to the closed universes contracting from the unstable de Sitter node towards the stable 
de Sitter node. 
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0.5 1 1.5 2 2.5 3 

X 

FIG. 3: Phase space for the system (|18(l reconstructed from the potential function (|16|) for the 
best fitted parameters (TABLE [HJ A3 = 0.3). The coloured domain of phase space is the domain 
of accelerated expansion of the universe. The red curve represents the flat model trajectory which 
separates the regions with negative and positive curvature. 



X 




FIG. 4: Phase space shown in Fig. |3] transformed on the compact Poincare sphere. Non-physical 
domain of phase space is marked as pink. 
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The trajectories situated in region III correspond to the models expanding towards stable 
de Sitter node at infinity. Similarly, the trajectories in the symmetric region II represent the 
universes contracting from the unstable node towards the singularity 

The main idea of the qualitative theory of differential equations is the following: instead 
of finding and analyzing an indiwidual solution of the model one investigates a space of 
all possible solutions. Any property (for example acceleration) is believed to be realistic if 
it can be attributed to a large subsets of models within the space of all solutions or if it 
possesses certain stability property which is shared also by all slightly perturbed models. 

IV. CONCLUSION 

The possible existence of the unknown form of matter called dark energy has usually been 
invoked as the simplest way to explain the recent observational data of SNIa. However, the 
effects arising from the new fundamental physics can also mimic the gravitational effects of 
dark energy through a modification of the Friedmann equation. 

We exploited the advantages of our method to discriminate among different dark energy 
models. With the independently determined density parameter of the Universe (fl m ,o — 0.3) 
we found that the current observational results require the cosmological constant n ~ in 
the Cardassian models. On Fig. Q we can see that both in the case of sample A (Fig. [Ifc) 
and sample C (Fig. ^i) n should be close to zero. Similarly if we assume that the density 
parameter for barionic matter is fl m ,o — 0.05 then n ~ 0.36 in the case of sample C (Fig. [I]p) 
and n ~ 0.39 for sample A (Fig. ^d). Moreover, we showed (for the sample C of Perlmutter 
SNIa data) that a simple Cardassian model as a candidate for dark energy is ruled out by 
our analysis if 0.36 < n < -0.50 for Q mfi = 0.3 and if 0.48 < n < 0.13 for Q mfi = 0.05 at 
the confidence level 2a. 

Therefore, the standard ACDM seems to be more appropriate for a large range of model 
parameters than the dark energy models based on the first version of Cardassian models. 

Our main result is that the structure of the phase space of accelerating models can be 
effectively reconstructed from SNIa data. 
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